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Abstract. Cross sections have been computed for rotational transitions of H2, 
' induced by collisions with H atoms, using a recent H - H2 potential calculated by 

Mielke et al. £P. These results are compared with those obtained with earlier 
potentials. Significant discrepancies are found with results deriving from the potential 
J> ■ of Boothroyd et al. J3] in the low collision energy regime. We compare also cross 

V) , sections derived using different levels of approximation to the vibrational motion. 
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The H - H 2 atom-molecule system is the simplest triatomic system; but, in spite 



of its structural simplicity, it has proved difficult to calculate the interaction potential 
to an accuracy which is sufficient to determine reliably the rotational excitation cross 
sections near threshold. In the regime of collision energies E <^ 5000 K, a potential 
barrier inhibits H-atom exchange scattering, and rotational excitation of the molecule 
occurs predominantly through non-exchange scattering; at higher energies, exchange- 
scattering occurs also. These processes contribute to the excitation of H 2 in the 
interstellar medium, notably in jets associated with star formation. Pure rotational 
transitions of H 2 have been observed from levels up to v = 0, J = 27 (Rosenthal et al. 



by means of the Infrared Space Observatory. 

A recent study of the H - H 2 potential by Mielke et al. [T] appears to have improved 
on the accuracy of earlier work by Boothroyd et al. |Sj and Partridge et al. [4J. Whilst 
all these calculations agree with regard to the general features of the potential, which 
becomes repulsive at short range and displays only a shallow Van der Waals minimum 
of approximately 2 meV (20 K), there remain significant differences in the predicted 
anisotropy of the potential; and it is the anisotropy which determines the magnitude 
of the rotationally inelastic cross sections. In this Letter, we present the results of 
calculations of cross sections for rotational transitions within the v = vibrational 
ground state of H 2 , comparing the values obtained using different determinations of the 
H - H 2 interaction potential. 

When determining the rotational excitation cross sections, it is customary to expand 

the interaction potential V(R,6,r), in terms of Legendre polynomials, P\(cos6), 

00 

V(R,9,r) = Y,vx(R,r)P x (cos9), (1) 

A=0 
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Figure 1. The isotropic vo and first anisotropic V2 potential expansion coefficients 
for the Boothroyd et al. and Mielke et al. potentials, computed in both cases for the 
equilibrium internuclear separation, r = 1.401 do, of the H2 molecule in its ground 
vibrational state. The unit of the potential is K; lE^ = 315780 K. 



where R is the vector connecting the centre of mass of the H 2 molecule to the H atom, 
r is the vector connecting the nuclei of the molecule, and cos 9 — R • f . As H2 is 
homonuclear, only those terms in the expansion with A even are non-zero. The results 
in Figure H were obtained using a least-squares minimization procedure, which provides 
an excellent fit to the original potentials when terms through A = 10 are included in the 
expansion. 

The behaviour of the term v 2 (R, r) in the classically allowed regions, where 
E ^ Vq (R, r), essentially determines the magnitude of the A J = 2 rotationally inelastic 
cross sections, in the range of collision energies considered here. Figure Q shows the 
isotropic term, Vq, and first non-zero anisotropic term, t> 2 , derived from the potential 
energies calculated by Mielke et al. (L and Boothroyd et al. jS] for the equilibrium 
internuclear distance, r = 1.401 ao, of the H2 molecule in its ground vibrational state, 
v = 0. There is excellent agreement for the isotropic term, v , but, as the expanded 
plot to the right in Figure Q shows, there are important differences for the first non-zero 
anisotropic term, t> 2 . The calculations of Boothroyd et al. imply a v 2 term which remains 
negative beyond its minimum, at R = 2.25 ao, whereas the calculations of Mielke et al. 
lead to a t>2 term which varies between positive and negative values for 3.7 < R < 12.9 
a . For atom-molecule separations, R, which correspond to low classically-allowed 
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Figure 2. Cross sections for the J = 2 — > transition, in units of 10 
centre-of-mass collision energy is expressed in K, relative to the J — level. The full 
curve was obtained with the potential of Mielke et al., and the broken curve with the 
potential of Boothroyd et al.: (a) at the equilibrium internuclear distance, r = 1.401 ao, 
in the vibrational ground state of the H2 molecule; (b) using the expectation values of 
V\ (i?, r) in the vibrational ground state (Equation ©). 



energies, the v-i term is small in magnitude for both potentials. Therefore, we expect the 
rotationally inelastic cross sections to be small at low collision energies. This expectation 
is confirmed by the results reported below. It is evidently a difficult task to determine 
reliably the value of the anisotropic term when its magnitude is as small as shown in the 
right hand panel of Figure ^ However, we note that, in the range of R which concerns 
us here, the Vi coefficient which derives from the calculations of Mielke et al. is in 
distinctly better agreement with the 1993 results of Partridge et al. [I] than with the 
1996 results of Boothroyd et al. [H|. 

As the collision energy increases and the distance of closest approach decreases, 
the magnitude of the i>2 coefficient increases and we find better agreement between the 
results deriving from the Mielke at el. and Boothroyd et al. potentials. Thus, we expect 
the rotationally inelastic cross sections to increase in magnitude, and the two potentials 
to yield similar results. Once again, we confirm this expectation below. 

In addition to calculating the cross sections with the internuclear separation fixed 
at its equilibrium value, r = 1.401 ao, in the vibrational ground state, we adopted a 



Letter to the Editor 



4 



more exact approach, consisting of evaluating 

Vx (R) = p ip* ( r ) v x (R, r) (r) dr (2) 

where ip(r)/r, the vibrational ground state eigenfunction at internuclear separation 
r, differs significantly from zero for r\ < r < r 2 . The vibrational eigenfunction was 
determined using the method of Marston and Balint-Kurti [Sj and the potential of 
Mielke et al. or Boothroyd et al., as appropriate. 

Rotationally inelastic cross sections were computed using the MOLCOL code [6 
and the exact coupled channels (CC) method 7\. In order to ensure convergence of the 
cross sections, for the range of collision energies considered in this Letter, we included 
the first 18 rotational states of the H 2 molecule, < J < 17, in the basis set. The J = 17 
level lies at an energy of 21412 K above J = 0. Because hydrogen-atom exchange is 
assumed not to occur, the ortho and para forms of H2 remain distinct, and separate 
calculations were performed including either the levels with J odd or the levels with J 
even. The first six terms (A = 0, 2, 4, 6, 8, 10) in the expansion of the potential were 
used in the calculations. 

In Figure EK, we compare the cross section for de-excitation, J = 2 — > 0, as 
calculated with the potentials of Mielke et al. [.tj and Boothroyd et al. [3], for the 
case of a fixed internuclear separation, r = 1.401 a . As anticipated, the cross-section 
tends to increase with the collision energy, and results obtained with the two potentials 
come into agreement. However as the collision energy decreases, the results obtained 
with the potential of Mielke et al. fall increasingly below those derived from the other 
potential, before rising again towards the J = 2 threshold. These differences relate to 
to the t>2 term in the potential expansion (see Figure Q). We have verified that the 
cross sections shown in Figure EK are reproducible by the independent scattering code 
MOLSCAT jH], to within the plotting accuracy. 

Figure EK shows the corresponding results for the J — 3 — > 1, transition of ortho- 
H 2 . The behaviour of this cross section is similar to that seen in Figure EK> for the 
equivalent transition in para-H2, except close to threshold (E = 1015 K), where both 
potentials yield a cross section which decreases with E. Once again, these results have 
been verified using MOLSCAT. 

In FiguresEb andEb, the rotational de-excitation cross sections for para- and ortho- 
H 2 are compared for the case in which Equation (j2J) was used to evaluate the potential 
expansion coefficients. Results obtained in this way should be more accurate than those 
derived assuming the internuclear distance, r, to be fixed at its equilibrium value. We 
see features in the cross sections which are similar to those apparent in Figures 
and EH- The discrepancies between the results obtained using different potentials are 
even greater than previously for collision energies E <; 3000 K. Additional calculations 
were performed for models in which simple harmonic oscillator and Morse oscillator 
eigenfunctions were used in Equation (J2J); they yielded results which are intermediate 
between those in Figures |2] and El and show the same features. 
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Figure 3. As for Figure but for the J = 3 — > 1 transition. The centre-of-mass 
collision energy is expressed in K, relative to the J = level. Cross sections obtained 
using v\(R, r = 1.401 do) are shown in panel (a), and those derived using Equation (J2J) 
in panel (b). 



As mentioned above, the behaviour with collision energy and the magnitudes of 
the cross sections relate to the t; 2 term in the expansion of the interaction potential. 
For scattering at low energies, the region R ^ 4 do is relevant, where the v% coefficient 
calculated by Mielke et al. [I] is in distinctly better agreement with that of Partridge et 
al. jlj than that of Boothroyd et al. j^j. Consequently, it was not surprising to find that 
the cross sections which derive from potentials pQ and jl] are also in better agreement. 
The rms and largest errors in the interaction energies computed by Boothroyd et 
al. (85.59 K and 1958.57 K, respectively) and Mielke et al. (3.25 K and 87.61 K, 
respectively) are to be compared with |u 2 | << 100 K for R ^ 4 a . Thus, the rms error 
quoted for surface of Mielke et al. [Tj approaches that required to calculate reliably the 
cross sections for at least the A J = 2 rotationally inelastic transitions. In the absence of 
any independent test of the accuracy of the calculated interaction potentials, it is clear 
that the cross sections computed with the potential of Mielke et al. are to be preferred. 

We note that the low-energy cross sections which derive from the potential surface 
of Mielke et al. [l] are even smaller than those obtained using the surface of Boothroyd 
et al., with values which fall to the order of 10~ 19 cm 2 for E « 2000 K; the corresponding 
thermal rate coefficient is of the order 10~ 13 cm 3 s _1 . This value is 2 orders of magnitude 
smaller than the rate coefficients for the J = 2 — > and J = 3 — > 1 transitions in H~2, 
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induced by collisions with He atoms at T = 2000 K j^j. In the interstellar medium, 
the elemental number density ratio of helium to hydrogen is nne/^H = 0.1. Thus, even 
if the hydrogen is mainly in atomic rather than molecular form, collisions with helium 
atoms will dominate the rotational excitation of the residual molecular hydrogen, owing 
to the much larger values of the He-H2 rate coefficients. 
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